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Kondo insulators are emerging as a simplified setting to study both magnetic and metal-to- 
insulator quantum phase transitions. Here, we study a half-filled Kondo lattice model defined on 
a magnetically frustrated Shastry-Sutherland geometry. We determine a “global” phase diagram 
that features a variety of zero-temperature phases; these include Kondo-destroyed antiferromagnetic 
and paramagnetic metallic phases in addition to the Kondo-insulator phase. Our result provides 
the theoretical basis for understanding how applying pressure to a Kondo insulator can close its 
hybridization gap, liberate the local-moment spins from the conduction electrons, and lead to a 
magnetically correlated metal. We also study the momentum distribution of the single-particle 
excitations in the Kondo insulating state, and illustrate how Fermi-surface-like features emerge as a 
precursor to the actual Fermi surfaces of the Kondo-destroyed metals. We discuss the implications 
of our results for Kondo insulators including SmBe. 

PACS numbers: 71.10.Hf, 71.27.+a, 75.20.Hr 


Quantum criticality in the vicinity of antiferromag¬ 
netic order is of interest to a variety of strongly corre¬ 
lated electron systems [l[. Heavy fermion systems oc¬ 
cupy a special place in this context [H, Q. These sys¬ 
tems are typically described by a Kondo lattice model, 
which contains a lattice of local-moment spins coupled 
to a band of conduction electrons. The strong interac¬ 
tions, which underlie the local moments, allow the ground 
states of these systems to be readily tuned by external 
parameters such as pressure or magnetic field. For this 
reason, antiferromagnetic (AF) quantum critical points 
(QCPs) have been identified in a host of heavy fermion 
metals. Experimental studies at such QCPs [j-lzl have 
provided evidence for a Kondo-destruction local QCP [H- 
ll0|, across which the Fermi surface jumps from “large” 
(incorporating the /-electrons) to “small” (excluding the 
/-electrons). Experimental efforts have also been devoted 
to heavy fermion metals that allow a systematic tuning of 
inherent quantum fluctuations, through strong magnetic 
frustration 15| . dimensionality [ibl or other related 


IThllOj. which are shedding light on novel phases 


that involve the development or destruction of not only 
the AE order parameters but also the Kondo entangle¬ 
ment. Such effects have been theoretically considered, 
and a global phase diagram has been advanced 

In the context of such developments in heavy-fermion 
metals it is then a natural question to ask, whether and 
how quantum phase transitions can be realized in Kondo 
insulators [2^? When the local-moment spins form a 
Kondo singlet with the spins of the conduction electrons. 


the resulting Kondo resonances are described in terms of 
a hybridization with the conduction electrons [25[. An 
incommensurate electron-filling implies that the chemi¬ 
cal potential intersects the hybridized bands, which leads 
to a heavy-fermion metallic state. By contrast, for a 
commensurate filling, the chemical potential falls in the 
middle of a hybridization gap m 113 , and a Kondo- 
insulator state ensues. If, by analogy with the case of 
the heavy fermion metals, an external parameter such 
as pressure tunes the system across a Kondo-destruction 
transition, the gap of the Kondo insulator will close. At 
the same time, the local-moment spins will be liberated 
from the conduction electrons, thereby yielding magnetic 
states in which the spin-rotational invariance is either 
spontaneously broken {e.g., an AE order) or preserved (a 
valence-bond solid or a spin liquid). This type of quali¬ 
tative considerations have led to a proposed global phase 
diagram for Kondo insulators [28|, but systematic theo¬ 
retical studies have yet to be performed to map this out. 

Studies along this direction are also important to un¬ 
derstand the on-going experiments on Kondo insula¬ 
tors [i^, which in recent years have been particularly 
fueled W the search for topological Kondo insulating 
states [29| in SmBe and related systems. Two types 
of experiments have added to the motivation for the 
present work. Eirst, it is known that applying pressure 
drives SmBe from a Kondo insulator to an AE metal , 
with indications for linear resistivity near the threshold 
pressure [3l|. Second, torque magnetometry measure¬ 
ments have observed a de Haas-van Alphen signal in 
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SmBe [32|, [33|, raising the exciting possibility that the 
Kondo insulator state harbors an incipient Fermi surface 
of the underlying conduction electrons (ssj. 

In this paper, we study a half-filled Kondo lattice 
model. First, in order to have a convenient tuning pa¬ 
rameter for quantum fluctuations of the local-moment 
magnetism, we focus on the model defined on the geomet¬ 
rically frustrated Shastry-Sutherland lattice. The calcu¬ 
lated phase diagram is strikingly similar to the global 
phase diagram proposed earlier based on qualitative con¬ 
siderations [HI. Second, we study the momentum dis¬ 
tribution of single-particle excitations in the Kondo in¬ 
sulating state. We identify a Fermi-surface-like feature 
in the Kondo insulating phase, and show how it rep¬ 
resents a precursor to the actual Fermi surfaces of the 
Kondo-destroyed metals that are “nearby” in the zero- 
temperature phase diagram. 

The half filled (z.e. one conduction electron per site) 
Kondo-Heisenberg model is defined as 


H = tij{c]„c 


ihj),cr 




-t-h.C.) + Jk JijSi ■ Sj 

i (ij) 

( 1 ) 

where (i^j) denotes a sum over neighboring bonds, with 
a hopping strength an RKKY interaction and a 
Kondo coupling The spin density of the conduction 
electrons at site i is = c\^{aai 3 / 2 )cii 3 ^ where CTais are 
the Pauli spin matrices, which are Kondo coupled (via 
Jk > 0) to spin-1/2 local moments, /2Yfis. 

To represent the spins, we use fermionic spinons [3^ fia 
which are subject to the constraint fjafio- = 1- The 
half filled condition is Uc = 1, where = ;^ Hi cr(4a^«cr) 
is the filling of the conduction band {N is the number of 
sites in the lattice). 

Global phase diagram: We consider the two- 

dimensional (2D) Kondo lattice model on the SSL ge¬ 
ometry, which is depicted in Fig. [TJ Here Ji and J 2 de¬ 
note the nearest neighbor (NN) and next nearest neigh¬ 
bor (NNN) RKKY interactions, respectively, and ti and 
t 2 are for the NN and NNN hoppings, respectively. In 
the limit of Jk = 0 both the Heisenberg model on the 
SSL [ 35 I, l36| and the electronic band structure have been 
studied [13, [HI in detail. It is known that for J 2 /J 1 > 2 
the SSL Heisenberg model exhibits an exact valence bond 
solid (VBS) ground state. In the incommensurate-filling 
case, quantum phase transitions between the heavy- 
fermion metal phase and other magnetic-metal phases 
have been studied [H, S, HH . 

We treat the Kondo interaction by introducing a 
Hubbard-Stratonovich decoupling in the Kondo-singlet 
{i.e. hybridization) channel Bi = Hcr^L/«cr* addi¬ 
tion, to capture the valence bond solid (VBS) and AF 
orders, we split the Heisenberg interaction into each of 
these respective channels written as Dij = Ha fLfjo- 
and = fj^{cr(^p/2)fig, and we weight the channels by 
a parameter x, where 0 < x < 1 [23|. This leads to the 
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FIG. 1. Zero temperature phase diagram as a function of 
frustration {J 2 /Ji) and Kondo coupling (Jx/ti), for half fill¬ 
ing ric = 1.0 and for 12 /ti = 0.25. The phases and their 
transitions are described in the main text. 


mean field Hamiltonian 

Hmf = C- Y.{b*clK + h.c.) + ^ Kflfia ( 2 ) 

i i 

(hj) 

where the sum over a is implied, and we have de¬ 
fined Jij = (1 — x)Jij. The constant term is C = 
H, (2|6,|VJk - AO+Hq,,o(2|Q^.f 
The hybridization is hi = JK^Bi)j2^ and the Hubbard- 
Stratonovich parameters in the resonating valence bond 
(RVB) singlet channel are Qij = xJij{Dij){2. Lastly, the 
AF order parameter is = (Si) [with the ordering wave 
vector Q = (tt, tt)]. We solve for these parameters as de¬ 
scribed in Ref. [23| using a four site unit cell labelled by 
X = with the corresponding four constraints 

and Kondo hybridization Ax and bx respectively. This 
leads to ten RVB parameters QxnQyiiQx+yi and Qx-v 
where i = 1 — 4 and the geometry is specified in Ref. [23| . 
We determine the zero temperature phase diagram for 
various different values of ^ 2 /^ 1 , and fix x = 0.7. This 
choice of x is guided by the magnetic phase diagram at 
Jk = 0, where taking x = 0.7 captures the quantum AF 
ground state (in the limit of J 2 = 0) [H]; taking this 
solution as a candidate state we have checked that our 
results discussed below are robust within a certain range 
of X. 

In Fig. [H we show the phase diagram for ^ 2/^1 =0.25, 
which reveals the three relevant phases. For small Jx/^i 
and J 2 /J 1 the model is in the AF phase, defined by 0 < 
|M I < 1/2. Here, Qij is only non-zero along the vertical 
and horizontal bonds and 6x = 0. As this phase has no 
Kondo screening and is antiferromagnetic we dub it AF s , 
where the S denotes a small Fermi surface. In the limit 
of large frustration and small Jk lA the model gives rise 
to the expected SSL-VBS, where Qx+y and Qx-y are the 
only non-zero singlet parameters. In the limit of large 
Kondo coupling we find a Kondo insulating (KI) phase. 
Here all the Qij and bx are non-zero, and they preserve 













3 



kjn 


FIG. 2. (Color online). Effects of the underlying Fermi sur¬ 
face on the band gap in the KI phase. The band gap AEk is 
plotted in momentum space for Jx/t — 1.5 in the KI phase. 
The minimum of the band gap is along the underlying small 
Fermi surface (the same as for Jk — 0.0), corresponding to 
the diamond marked in red. 


the symmetry of the Shastry-Sutherland lattice. 

We find direct transitions between AF s and KI, as well 
as between SSL-VBS and KI. In the incommensurate- 
filling case [m, we found various intermediate phases 
(between the VBS and heavy Fermi liquid) that break the 
underlying symmetry of the Shastry-Sutherland lattice 
and as a result exhibit partial Kondo screening. Here, 
we find these solutions to be energetically not compet¬ 
itive. While our approach yields first order transitions 
(see the Supplementary Material [4l|), it is important to 
see how the fluctuation effects beyond our approach affect 
the nature of the transitions. Based on the studies [^| 
on the quantum fiuctuations in pertinent quantum impu¬ 
rity models in the context of extended dynamical mean 
field approach 311^, we expect the transitions to be 
continuous. 


Two remarks are in order. First, the phase diagram 
we have found is strikingly similar to the global phase 
diagram proposed for Kondo insulators on qualitative 
grounds [28|. From a more general point of view, we can 
think of J 2 / Ji as a measure of the strength of quantum 
fiuctuations in the system. Therefore, the phase diagram 
we have derived is representative of Kondo insulating sys¬ 
tems more generally; for instance, dimensionality tun¬ 
ing could serve for a similar purpose [l^. Second, each 
type of phase transition between the KI phase and either 
the AF^- or the SSL-VBS is actually a metal to insula¬ 
tor transition. Thus, there should be significant effects 
on Fermi surface probes as the transition is approached 
from the metallic side as well as from the insulating side, 
on which we now turn to. 


Single-particle excitations and precursor of the small 
Fermi surface in the KI phase: We focus on the KI 
phase itself to consider the momentum distribution of 
the single-particle excitations. We are primarily inter¬ 
ested in what effects the underlying Fermi surface of the 
conduction electrons have “imprinted” on the properties 



FIG. 3. (Golor online). The spectral function for Jx/ti = 1.5 
along the content momentum cuts (a) k — ky — kx and 
(b) /c = /cx, /cy = 0, the red lines mark the location of the 
small Fermi momentum kp- Energy dispersion curves of the 
spectral function are shown along the momentum cuts (c) 
k — ky — kx and (d) k — kx^ ky — t). Here, each curve is 
shifted vertically by an integer corresponding to the wavenum¬ 
ber index and we have broadened the delta function by a 
Gaussian of width O.lt. 


of the insulating state. For definiteness, we will denote 
such small Fermi momenta by kj?, and the large Fermi 
momenta (those that incorporate the /-electrons, see be¬ 
low) by k.p. To make concrete connections, we now re¬ 
strict ourselves to a model defined on a 2D square lattice 
with NN couplings only, and therefore consider one site 
per unit cell. In addition, for the KI phase where the 
magnetic order parameter vanishes, we no longer con¬ 
sider the effects of the Heisenberg term in Eq. m and 
set Ji = J2 = 0; this still keeps the salient properties of 
the momentum distribution in this phase. Note that in 
the following we only have one hopping parameter t and 
we drop the subscript “1” in the remainder of the paper. 

We first consider the hybridization gap AE\^ as a func¬ 
tion of momentum. Solving for the band structure in the 
square lattice case yields two bands E]^± = ^(ck — A) ± 

V where Ck = —2t{coskx + cosky) — ja is 

the dispersion with a bandwidth W = 8t and chemical 
potential fi. A Kondo insulator of course has no Fermi 
surface. However, plotting the direct hybridization gap 
AF^k = F^k+ — F^k- as a function of k in the entire Bril- 
louin zone reveals a special surface (line) in the Brillouin 
zone. As shown in Fig. [21 AF^k is minimized along the 
diamond marked in red, which corresponds to the small 
Fermi surface for the half-filled conduction-electron band, 
TT-c = 1. We therefore reach one of our main results, 
namely AF^k is minimized on the small Fermi momenta, 
ki?. The magnitude of the direct (and indirect) gap is 
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FIG. 4. (Color online). Momentum distribution of the c elec¬ 
trons rifc, for k — kx — ky. rik versus k for various Kondo 
couplings in the KI phase (a) and the metallic phase with 
Tic = 1.25 (b). The overall shape across /cf in both cases is 
discussed in the text, while for the metal there is a jump at 
the large Fermi wave-vector kp (dashed line) of the size of the 
quasiparticle residue. 


discussed in the Supplementary Material (Ref. 0). 

To further illustrate the role played by the small 
Fermi momenta ki?, we consider the dispersion of the 
single particle excitations and momentum evolution of 
the conduction-electron spectral function = 

—ImGc(k, co’)/7r (see the Supplemental Material [^) as 
shown in Fig. [3l Along the cut kx = ky^ as shown in 
Fig. [3](a), we find the quasiparticle states dispersing to¬ 
wards the small Fermi momentum k^ = (7r/2,7r/2) as 
the Fermi energy is approached, although they are even¬ 
tually gapped out by the hybridization when k gets too 
close to ki?. This point is also illustrated in the energy 
dispersion curves, shown in Fig. [3](c). The same trend 
is also observed for the momentum cut ky = 0, Figs. [3] 
(b) and (d), although here the small Fermi momentum is 
ki? = (tt, 0), which is located on the zone boundary. 

To appreciate the above observations, we note that 
the small Fermi momenta k^? are special because, for 
Jr = 0, the momentum distribution of the c electrons, 
rik = cr^k,cr )5 hus u jump of exactly 1 across such 

momenta. A non-zero Kondo coupling will smear this 
jump [m, but this smearing occurs gradually. Indeed, 
as shown in Fig. [U (a), near the small Fermi momentum 
ki?, for various values of the Kondo coupling in the KI 
phase, we see a step function for Jk = 0 develop into 
an “S-shape” pinned at ki? with increasing Jk- (With¬ 
out a loss of generality, we consider the trajectory in the 
Brillioun zone 0 < kx = ky < ir.) This smeared jump in 
the momentum distribution of the occupation number at 
ki? is caused by the same physics that induces the small 
excitation gap at ki? illustrated in Figs. [2] and [3l 

The simplicity of the Kondo insulator is that the 
large Fermi momenta, k^, are located at the Brillouin 
zone boundary. This is to be contrasted with the 
incommensurate-filling case, as illustrated in Fig.lU^b) for 
TT-c = 1.25. In this case, the features at the small Fermi 
momenta, k^, remain similar to the Kondo-insulator 
case. However, now, k.p occurs in the middle of the 
Brillouin zone. As is characteristics of a heavy-fermion 
metal [25[, nk display a sharp drop at k.^. The drop 



Pressure 


FIG. 5. (Golor online). Schematic finite temperature phase 
diagram as a function of the control parameter 6. The close 
proximity of the Kondo-insulator (KI) state to a Kondo- 
destroyed AF phase (AF^) gives rise to a zero temperature 
QGP and a concomitant quantum critical (QG) regime. 


is tiny, measuring the exponentially small quasiparticle 
residue. 

Discussion and outlook: We now explore the impli¬ 
cations of the global phase diagram shown in Fig. [H By 
tuning various Kondo insulating compounds (e.g. un¬ 
der pressure) into and out of various ground states, this 
phase diagram opens up the exciting possibility of re¬ 
alizing new types of quantum phase transitions. From 
the Kondo-insulating state, such tuning can suppress the 
insulating gap by destroying the Kondo effect, liberate 
the local-moment spins, and give rise to either an AF or 
paramagnetic metallic phase. In these metallic phases, 
the Fermi surface is small as defined earlier, namely it 
incorporates only the conduction electrons but not the 
/-electrons. 

Taking a cut in the global phase diagram provides the 
basis to construct a schematic finite temperature phase 
diagram. We illustrate this point in Fig. [5] on the case 
when, along the chosen cut, the Kondo-destroyed phase 
nearby to the Kondo insulator is AF^. The AF order 
is destroyed at the Neel temperature Tw, whereas the 
Kondo insulating state persists up to a renormalized ef¬ 
fective Kondo energy scale T*, which is distinct from the 
much higher in energy bare single ion Kondo tempera¬ 
ture T^, due to a renormalization from the presence of 
the RKKY interaction. Thus at finite temperature, by 
tuning an external control parameter ((5) the system can 
go from an AF s phase into a cross over regime driven by 
the zero temperature quantum critical (QC) fluctuations, 
which then crosses over to the KI phase at T*. Starting 
from high temperature for a control parameter that is in 
the KI phase at zero temperature, lowering the tempera¬ 
ture through leads into the QC regime and non-trivial 
scaling features, which then become cut off by the Kondo- 
insulator state below T*. This provides the theoretical 
basis to understand why pressure induces a Kondo insu- 
lator to antiferromagnetic metal transition 


31 ,^ 


It also makes it at least plausible that a precursor Fermi 
surface appears in a Kondo insulator, similar to an inter- 
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pretation of the dHvA effect measured in SmBe [^ . 

In conclusion, we have determined the global phase di¬ 
agram in the Kondo-insulator regime of a Kondo-lattice 
model. This phase diagram makes it natural that pres¬ 
sure tuning of Kondo insulators induces magnetic metal 
phases. Our results should also motivate the study of an¬ 
tiferromagnetic and insulator-metal quantum phase tran¬ 
sitions in Kondo insulators with varying geometrical frus¬ 
tration, or with varying dimensionality through thin films 
or heterostructures. Finally, we have studied the momen¬ 
tum distribution of the single-electron excitations in the 
Kondo-insulator phase, and demonstrated the imprints 
of a small Fermi surface in this distribution. 
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SUPPLEMENTARY MATERIAL 

Here we discuss the various equations we have used in studying the nature of the Kondo insulating phase. We show 
the mean field parameters as a function of the Kondo coupling Jk across each phase transition in Fig. [HD (as in the 
phase diagram Fig. 1 of the main text) and also discuss the direct and indirect single-particle excitation gaps in the 
Kondo insulating phase. 


EQUATIONS FOR THE SINGLE PARTICLE EXCITATIONS IN THE KONDO INSULATING PHASE 


Here we give the equations for the occupation number and spectral function within our approach. In the absence of 
the RKKY interaction focusing on a square lattice, and after the Hubbard-Stranotovich decoupling, the Hamiltonian 
becomes 

Hmf = X (/I/- - ^) + E (4aC,> + h.c.) + A ^ ^ {hlcln, + h.c.) . 

Focusing on translationally invariant solutions, we can set iXi = —A and bi = b. We can then diagonalize the 
Hamiltonian via Fourier transform and obtain 


Hmf = He + ^ ^^k+7ka+7ka+ + ^k-7ka-7/ca- 
k,cr 


(SI) 


where the light and heavy bands (+/—) are 


= l(ek - A) ± 

and 

The quasiparticle operators, 7kcr±: are related to the original degrees of freedom via 


Ckcr — '^k 7 kcr-|- T '^k 7 kcr— 5 

/kcr ~ '^k 7 kcr-|- '^k 7 kcr— 5 


where -h = 1 and the form factors are 


^^2 _ 1 ^ 1_ ^k + A _ 

" 2 + 2y(ek + A)2 +4P’ 

^ 2 _ 1 1 Ck + A 

2 ■ 27 (ek + A)2+462- 


With these results in hand we can now determine the Green function of the conduction electrons, 


^c(k, ioJn) — . 


lUJrn 


■E, 


'k+ 


lUJrn 


^k- 


(S2) 


(S3) 


(54) 

(55) 


(56) 

(57) 


(S8) 
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with a spectral function 


Ac(k, w) = ul^5{uj - Eu.+ ) + vl,5{uj - E’k-)- (S9) 

Consequently, the occupation number of the c electrons at zero temperature is 

nk = ulQ{-E^,+) + ^;^0(-Ek-), (SIO) 

where 0(x) is the Heaviside step function. From the Green function we can also determine the quasiparticle residue 
in the heavy Fermi liquid phase for the filling ric = 1.25 as in the main text. In the heavy Fermi liquid phase the total 
filling will he Tic rif = 2.25, hence = 0 and there is only a pole at away from zero energy. Therefore, 

the residue for the c and / electrons at the Fermi energy is and Zf = respectively. 


MEAN FIELD PARAMETERS 


(a). 

-J 2 /J 1=0.40 
t2/ti=0.25 

Mq -.-.T™- 
Qxi/Ji .■. 

QxV^i . • . ■ 

byjK 





• • 



0.6 0.8 1 1.2 1.4 

Jk/Ii 


0.6 

(b) 

' Mq' .-t:. 

0.5 


• Q,,i/Ji .■. 

0.4 

J 2 /Ji= 1.60 

Qx+y/Ji . • . 

1 byjK 

0.3 

t 2 /ti= 0.25 


0.2 



0.1 




0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 

Jk/Ii 


FIG. SI. Mean field parameters across each of the Kondo insulating to Kondo destroyed transitions for t 2 = 0.25ti. (Only the 
distinct parameters are shown.) (a) For the AFs to KI phase transition as a function of the Kondo coupling Jk for J 2 I Ji — 0.4. 
We find a direct first order transition between the AF s and KI phase signaled by a vanishing of the magnetic order parameter 
Mq (where Q corresponds to the AF wave vector); (b) Mean field parameters for the SSL-VBS to KI phase transition as a 
function of the Kondo coupling Jk for J 2 /J 1 = 1.6. We find a direct first order transition between the VBS and KI phase. 


DIRECT AND INDIRECT GAPS 

We define the direct gap (AE* = miuAEk), which captures the lowest lying excitationsthat don’t change momentum. 
Likewise, we define the indirect gap {dE = minEk+ — maxEk-), which corresponds to the lowest-energy excitation 
when the momentum is allowed to change. In Fig.[S2](a), we show the direct gap as a function of the Kondo coupling 
Jx, which has the expected exponential dependence. Compared to the indirect gap, which is of the order of the 
Kondo temperature, the direct gap is larger and has a square-root dependence on the latter, as is shown in Fig. [S2] 
(b). 
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FIG. S2. Direct and indirect gaps in the Kondo insulating phase, (a) Direct gap AE as a function of the Kondo coupling Jk- 
Here, the dashed line is a fit to the expected exponential form [aexp(—5 /Jk)]. (b) The direct gap as a function of the indirect 
gap 6E, satisfying AE ~ VSE. 




